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Abstract 



K.^ Exact diagonalization studies for a doped t-J ladder (or double chain) show 

o 

\f^ • hole pairing in the ground state. The excitation spectrum separates into a 



limited number of quasiparticles which carry charge +|e| and spin ^ ^^id a 
triplet mode. At half-filling the former vanish but the latter evolves continu- 



C^ ' ously into the triplet band of the spin liquid. At low doping the quasiparticles 

a" 

form a dilute Fermi gas with a strong attraction but simultaneously the Fermi 



wavevector, as would be measured in photoemission, is large. 
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The properties of strongly correlated electrons confined to a ladder (or double chain) 
and described by t-J or Hubbard models have been the subject of intensive investigation 
recently l|l]-§]. The reason lies in the unusual spin liquid nature of the undoped parent 
system |1],0J^. The key question is the evolution of the finite gap in the spin excitation 
spectrum upon doping. The spin gap remains in other spin liquids systems and is a sign 
of strong superconducting fluctuations i9|JT0(|. A recent analysis of the t-J ladder using a 



mean field (MF) theory with Gutzwiller renormalization of the matrix elements to account 
for the strong correlations, gave a continuous evolution of the spin gap with doping [^ . The 
short range resonance valence bond (RVB) state evolves into a superconductor with modified 
(i-wave symmetry within this MF theory for the t-J ladder. A tendency towards modified 
(i-wave superconductivity was also found in a renormalization group calculation and in a 
very recent numerical study of the Hubbard ladder although no actual enhancement of the 
superconducting correlations was reported Q. In this letter we report on the properties of 
finite t-J ladders up to a size of 10 x 2 sites using a Lanczos diagonalization method. We 
find results which give clear evidence of hole pairing and a modified (i-wave RVB state in 
lightly doped systems in agreement with the MF theory. An interesting difference however 
is the discontinuous evolution of the excitation spectrum upon doping. Upon doping new 
quasiparticle (QP) excitations appear carrying both charge and spin. These excitations are 
in addition to a band of spin triplets which evolve continuously away from the undoped 
spin liquid. This separation of the excitation spectrum into bound holon-spinon QP's and 
collective triplet excitation contrasts with the full spin-charge separation found in a Luttinger 
liquid. 

Another reason for especial interest is the possibility of realizing a lattice of weakly 
coupled ladders in the compounds SrCu20-i [|TT| and also (V"(9)2-P2C'7 |T^. In the first case 
the ladders are isotropic with equal hopping matrix elements and exchange couplings along 
rungs and legs and we concentrate on this case. 

The t-J ladder we will study reads 



^ = -^EiE cL(j)Ca.(j + l) + cL(j>2.(j)] + H.C.} 
j,cr a=l 



J,a 



+J'E[Si(j)-S2(j)-^i(j>2(j)], (1) 

i 

where j runs over L rungs, and a (=1, J.) and a (= 1, 2) are spin and leg indices. We take 

t = 1 as the energy units. The first term is the kinetic energy and the J (J') are exchange 

couphngs along the ladder (rungs). The local constraint excludes double occupancy at every 

site {nai;{j)nai{j) = 0). Although the isotropic case J' = J is of most interest, we study also 

the limit of J' ^ J, 1, which can be easily understood and is a good starting point to trace 

back to the isotropic case. Periodic or antiperiodic boundary conditions (PBC, APBC) are 

used along the ladder and then the wavevector k = {k^, ky) is well defined, k^ (ky) are along 

the ladder (rungs). 

To begin we summarize the properties of spin excitations at half-filling. In this case (|I]) 
reduces to the Heisenberg spin ladder (aside from a constant). Several different approaches 
have shown it has a spin liquid ground state with a finite excitation gap |l|J^,P]. In the limit 
of J' — i> oo, the ground state is an ensemble of singlet rungs and the total spin is trivially 
S = 0. Low lying excitations are Bloch states of one triplet with energy J' + Jcosk^, and 
have a minimum at k^ = vr. With decreasing J'/ J, the dispersion becomes more linear with 
a smaller gap, but at J' = J the gap remains finite (Agpin = 0.5 J). 

The first evidence of enhanced superconducting fluctuations in the t-J ladder is the 
formation of a bound pair of two holes doped into half-filling. The ground state wave 
function is calculated by the Lanczos method up to L = 10 and found to have S* = and 
k = (0,0). The two holes have a positive binding energy defined by, E-q = 2Egs{'^L—1) — 
Eqq,{2L) — Egs{2L—2) , where Egs{N) denotes the ground state energy for A^ electrons. In 
Fig. |l], we show E-q for L = 8 and J = 0.3 as a function of the rung exchange, J'. The 
binding energy is positive down to the isotropic value, J' = J = 0.3 appropriate for a 
cuprate. The same calculation for four holes indicates no phase separation up to the largest 



value J' = 3.0. 

In the limit of strong rung coupling J' ^ J,l, this binding can be easily understood. At 
half-filling every rung forms a spin singlet. A single hole breaks a singlet so the second hole 
is attracted to the same rung to avoid breaking another singlet bond. The binding energy 
is therefore Eb ~ J' — 2 asymptotically. 

The hole binding in the ground state is also seen directly in the hole-hole correlation 
function, {nah{j)na'h{j')) i^ahU) = 1 — na{j): hole density). We calculated the size of 
bound hole pair ^ by fitting {nih{i)n2h{j)) ~ const, x (e"'*^-''/^ + e^^^^'*~-''-'/^) for \i—j\ = ■§ 
and ■§ — 1- The results are shown in the inset of Fig. [l| for L = 6, 8, and 10. In the strong 
coupling limit J' ^ J,l, the two holes are confined to the same rung and ^ ^ 0. The size 
of the bound hole pair increases with decreasing J', but is still quite short, ^ ~ 2, even for 
isotropic coupling J' = 0.3. Note at this value, the hole- hole correlation is maximum for 
holes on neighboring rungs, i.e., {nih{j)n2h{j + !))• 

One of the most interesting properties of the t-J ladder is that there are two distinct 
types of spin excitations upon doping. Let us start from the large J' limit, where they are 
easily distinguished. The first is the triplet excitation similar to that at half-filling and arises 
when a singlet rung away from the hole pair is excited to a triplet. The triplet propagates 
along the ladder with the matrix element -I, and can also pass through the hole pair. For 
the second type of spin excitation, the presence of the holes is essential. The bound hole 
pair in the ground state dissociates into two separate holons, each of which is now bound 
to a spinon on the same rung to form a QP with charge +|e| and spin | [|13[. In this sense. 



these holon-spinon bound pairs are similar to QP's in conventional Fermi liquids. 

The energies and allowed numbers of the two types of spin excitations are different. The 
triplet can be excited only at the rungs without holes, whereas the second type needs a hole. 
Therefore the number of possible excitations is proportional to (1 — 6) and 6, respectively, 
where 6 is the hole doping. 

The second type of spin excitations have lower energies in the large J' region as follows. 
For both types at least two singlets have to be broken with energy cost J', but their kinetic 



energies differ. In the first type, rungs are occupied by a hole pair and by a triplet. The 
hole pair and the triplet move independently along the ladder with the matrix elements 
(J' — jt)~^ and Y, respectively. In the second type, two rungs are occupied by only one 
electron in a bonding orbital with kinetic energy gain (—2). These singly occupied rungs 
move with the matrix element |. Therefore the total gain of kinetic energy is much larger 
for the second type when J' ^ 1, J. Accordingly the lowest spin triplet excitation is of the 
second type as shown in Fig. |l|, and Aspin ~ Eb- With decreasing J', the two excitations 
become closer to each other in energy and they are strongly mixed since they can have the 
same symmetry. 

The different characters of the two types of excitations are clearly distinguished by the 
hole-hole and spin-hole correlation functions. Figure ^ shows the correlation functions cal- 
culated by the Lanczos method for the two excitations with S"^ = 1 for J = 0.3 and J' = 0.3, 
and 3.0. In the QP excitation (dotted line), the hole-hole correlation is maximum for the 
longest separation, but for the spin-hole correlations on the same rung. This tendency is 
obvious for the large rung coupling J' = 3.0, but still holds down to the isotropic value 
J' = 0.3. For the other excitations (dashed line), at J' = 3.0 the two holes are strongly 
bound on the same rung and repel the triplet. With decreasing J', the hole pair becomes 
more extended as in the singlet ground state. The most dominant configuration at J' = 0.3 
is that the two holes are in the neighboring rungs and different legs. 

The two types of excitations have different contributions to the spin susceptibility and 
structure factor, 

Im Sik,uj) = J2\{a\S^\GS)\' 5iuj -E^ + Egs). (2) 

a 

Here |a) and |GS) denote an 5 = 1 eigenstate and the ground state with energies Ea and 
Eqs- As was discussed before the spin excitations with separate holes have lower energies. 
However at 5 ^ 1, their number is smaller (oc 6) than the excitations with the bound hole 
pairs and triplet rungs (oc 1 — 6). Therefore with decreasing temperature the susceptibility 
will show a large exponential drop at T ~ J' and a smaller drop further at T ~ Aspin, 



corresponding to the different excitation energies. Figure |^ shows the structure factor, 
Im S(k,uj), for two holes in 8 x 2 sites. Large peaks are seen around k = (vr, tt) and the 
energy of the order of J'. These are due to the excitations away from bound hole pairs. 
On the other hand, the QP excitations do not have large weights in spite of their lower 
excitation energies, because of the change in the charge configuration from the ground state. 
This point would be important when comparing to neutron scattering experiments. 

Finally we discuss the one-particle Green's function where we can see QP excitations 
directly. By using the Lanczos method combined with a continued fraction method, we 
calculate the spectral function for two holes. The electron and hole parts of the spectral 
function are defined as 

/le(k,cu)=5:i(a,2L-l|cL|GS,2L-2)|2 

a 

x6{u - E^{2L-l) + ^Gs(2L-2) + /i), 
A(k, uj)= Y, I (a, 2L-3|ck.|GS, 2L-2) p 

a 

x6{uj + E^{2L-3)-EGs{2L-2)+fi), (3) 

where \(y,N) is an eigenstate for A^ electrons with the energy Ea{N) and GS denotes the 
ground state. Positive (negative) energies correspond to the electron (hole) part. The 
chemical potential is defined by /i = |[£'gs (2-^—1) — -E'gs(2-^— 3)]. The results are shown in 
Fig. ^ for L = 8 and J = J' = 0.3. The wavevector along the ladder k^ = ^n {n: integer) 
are for PBC and k^ = "^{n + \) for APBC. The ground state energy EgsC^L - 2) and the 
chemical potential fi in Eq. (|]) are the average over both boundary conditions. 

There are large weights for the bonding {B){ky = 0) and antibonding {A){ky = it) orbitals 
only near the Fermi energy u = 0, and they seem to constitute QP bands. Away from the 
Fermi energy, the individual QP peaks are much less prominent and there is an incoherent 
part with the energy of the order of 1. 

The QP spectrum is consistent with the MF theory based on the rf-wave RVB state. 
Without holes the B- and A-bands are degenerate because the kinetic energy is zero. Upon 
hole doping, {ci^ij)c2aij)) on a rung becomes finite and the two bands are split. The B- 

6 



band is pushed down and occupied by more electrons, while the A-band is pushed up and 
occupied by less electrons. The QP with energy closest to uj = has a wavevector nearest 
to the original Fermi kp'- {kx = ^ for B and k^ = y' ^^ ^)- Because of the band splitting, 
k^ > kp, but the Luttinger sum rule is satisfied, k^ + kp = {1 — 6)11. This means the Fermi 
surface is large, consistent with photoemission experiments on cuprates. 

It is important to notice the QP peaks near the Fermi energy have their counterparts 
on the opposite side of the Fermi energy. An electronic QP peak at energy u > has a 
shadow hole peak at energy around —u < 0, and vice versa. These peaks indicate that 
the QP excitations are those of the Bogoliubov QP's as in BCS theory, i.e., mixture of an 
electron and a hole (aj^ = WkCk^ + fkC-k;)- The weights in the electron and hole parts are 
proportional to l^kP and |fkP- They are hole-like around k^ = and electron-like around 
kx = IT for both B and A bands. There exists a finite energy gap in the QP spectra. The 
electron and hole branches both come close to the Fermi energy at fc^; ~ ^, but instead of 
passing through they move away from it. The energy gap is 0.18t at k = (^, 0) and 0.23t at 
k = (^, tt), corresponding to a QP gap 2Aqp ~ 0.2(~ ^). It is interesting to note that the 
Stephan-Horsch results for A(k, cj) in two-dimensional clusters |T^ show similar behavior 
for k-points not along (1,1) but no shadow peaks for k || (1,1), indicating (i2:2_j^2 -pairing 
also. 

The results of our exact diagonalization studies show that the doped t-J ladder belongs 
to a different universality class than the single chain. Here the separation is into a limited 
number of quasiparticle excitations carrying charge +|e| and spin 1 which vanish as 6 —^ 0, 
and a triplet collective mode which remains as 5 — i> 0. So the triplet mode is not a collective 
mode of the quasiparticles. This separation resembles in certain ways the proposal of Sokol 
and Pines for a doped quantum critical regime [|T5| and also Chubukov and Sachdev for 
collective and Fermion contributions to the susceptibility |T6|. Note the quasiparticles do 
not behave as a usual Fermi liquid. Although their number is strictly limited when 6 <^ 1, 
yet their pairing energy is finite as 5 — > 0. In this sense they resemble a dilute Fermi gas 
with a strong attraction giving pair binding |l^. However if we look at the dispersion energy 



of an added hole that would be seen in photoemission, then from Fig. ^ it is characterized 
by the large Fermi surface. The quasiparticles do not form a usual Fermi liquid but display 
a new and interesting mixture of dilute Fermi gas and large Fermi surface behavior. 

The authors thank S. Gopalan, H. Monien, D. Poilblanc, D. Wiirtz, F. C Zhang, and 
especially M. Sigrist for fruitful discussions. The work was supported by the Swiss National 
Science Foundation Grant No. NFP-304030-032833 and SNF-21-27894.89 and by an internal 
grant of ETH-Ziirich. 



REFERENCES 

[1] E. Dagotto, J. Riera, and D. J. Scalapino, Phys. Rev. B 45, 5744 (1992); T. Barnes et 
al., Phys. Rev. B 47, 3196 (1993). 

[2] T. M. Rice, S. Gopalan, and M. Sigrist, Europhys. Lett. 23, 445 (1993). 

[3] M. Sigrist, T. M. Rice, and F. C. Zhang, preprint, ETH-TH/93-35. 

[4] R. M. Noack, S. R. White, and D. J. Scalapino, preprint. 

[5] S. Gopalan, T. M. Rice, and M. Sigrist, preprint, ETH-TH/93-39. 

[6] D. V. Khveshchenko and T. M. Rice, preprint; D. V. Khveshchenko, preprint. 

[7] R. Hirsch, Diplomarbeit, University of Koln (1988). 

[8] S. P. Strong and A. J. Milhs, Phys. Rev. Lett. 69, 2419 (1992). 

[9] M. Ogata, M. U. Luchini, and T. M. Rice, Phys. Rev. B 44, 12083 (1991). 

[10] M. Imada, Phys. Rev. B 48, 550 (1993), and references therein. 

[11] M. Takano et al., J JAP Series 7, 3 (1992). 

[12] D. C. Johnston et al., Phys. Rev. B 35, 219 (1987). 

[13] Similar excitations are also found for a t-J chain with frustrated couplings. K. Takano 
and K. Sano, Phys. Rev. B 48, 9831 (1993); I. Bose and S. Gayen, Phys. Rev. B 48 
10653 (1993). 

[14] W. Stephan and P. Horsch, Phys. Rev. Lett. 66, 2258 (1991). 

[15] A. Sokol and D. Pines, Phys. Rev. Lett. 71, 2813 (1993). 

[16] A. V. Chubukov and S. Sachdev, Phys. Rev. Lett. 71, 169 (1993). 

[17] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, Phys. Rev. Lett. 71, 3202 
(1993). 



FIGURES 
FIG. 1. Binding energy of two holes, spin gap, and energy of the triplet excitation away from 

the bound hole pair. Inset shows the size of the bound hole pair for the ground state 

FIG. 2. Hole-hole ((a),(c)) and spin-hole {{b),{d)) correlation functions for the two triplet 
states. Dotted lines are for the lowest triplet state (QP excitation) and dashed lines are for the 
triplet state with a bound pair. 

FIG. 3. Spin structure factor, Im 5(k, to), for the t-J ladder with two holes and APBC. L = 8 
and J = J' = 0.3. The width of each line represents the strength of the excitation. 

FIG. 4. Spectral function of the one-particle Green's function, A{\i.,uj), for two holes. L = 8 
and J = J' = 0.3. The width of each line represents the strength of the excitation. 
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